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AERODYNAMICS  OF  CONES  AT  HIGH  ANGLES  OF  ATTACK  IN  LOW  SUB- 
SONIC  FLOW 


SUMMARY 

This  report  is  a  theoretical  study  of  various  phenomena  observed  in  the  low  subsonic 
regime  when  a  slender  cone  of  revolution  is  pitched  to  high  angles  of  attack.  Three- 
dimensional  boundary  layers  are  calculated  up  to  the  line  of  primary  separation.  The 
structure  of  the  separated  wake  is  obtained  with  a  two-dimensional  model  and  the 
results  are  used  as  input  in  a  leeside  boundary  layer  calculation.  This  calculation  is 
carried  on  to  the  line  of  secondary  separation.  The  positions  of  separation  are  com¬ 
pared  with  experimental  results.  The  wake  is  found  to  be  unstable  at  the  angles  of 
attack  where  strong  asymmetric  forces  are  observed.  A  stability  analysis  of  the  boun¬ 
dary  layer  rules  out  streamwise -directed  vortices  as  a  contributor  to  the  asymmetric 
forces. 


1  INTRODUCTION 


Figure  1.1  Vortex  formation  on  the 
leeside  of  inclined  body 
of  revolution 

a)  attached,  b)  separated  symmetric 
and  c)  steady  asymmetric  flow 


Allan  and  Perkins  (1)  have  classified  flow 
around  a  slender  body  in  four  distinct 
groups  depending  on  the  angle  of  attack 
(Figure  1.1).  At  small  angles  the  flow  is 
attached.  When  the  angle  is  increased,  the 
boundary  layer  separates  on  the  leeside, 
producing  a  symmetric  vortex  system.  If  the 
body  is  pitched  to  even  higher  angles,  the 
symmetric  structure  breaks  down.  The  vor¬ 
tex  system  is  then  steady,  asymmetric  and 
produces  often  undesired  lateral  forces  even 
at  zero  sideslip.  At  very  high  angles  of  attack 
the  wake  becomes  unsteady. 

These  phenomena  are  well  known,  but  not 
well  understood.  The  goal  of  the  present 
thesis  is  to  shed  light  on  the  flow  patterns 
mentioned,  with  special  attention  to  steady 
asymmetric  flow.  The  unsteady  phase  will  be 
just  superficially  treated. 

Within  the  aircraft  and  missile  industry, 
where  the  present  requirements  for  ma¬ 
noeuvrability  are  difficult  to  meet,  substan¬ 
tial  research  is  carried  out  on  slender  bodies 
at  high  angles  of  attack.  A  number  of 
experiments  reveals  connections  between  the 
asymmetric  forces  and  important  flow  and 
geometry  parameters.  We  shall  now  review 
the  most  important  parameter  effects. 


1.1  Geometrical  imperfections 


Experiments  with  varying  surface  finish  have  been  performed  on  a  tangent  ogive 
forebody  (2).  The  surface  roughness  was  found  to  have  a  large  effect  on  vortex  induced 
side-force. 

When  a  body  is  tested  in  various  roll-positions,  the  side-force  may  change  direction  for  a 
range  of  roll  angles  (2,3, 4, 5).  There  is,  as  expected,  no  mirror  symmetry  in  the  results. 
Similar  results  are  obtained  when  only  a  removable  nose  is  rolled  (5).  This  strongly 
indicates  that  nose-tip  imperfections  are  the  triggering  mechanism  for  the  asymmetric 
forces. 


1.2  Mach  number 

From  various  experimental  works  (2,3,5),  the  following  general  trend  can  be  outlined. 

The  side-force  is  most  dominant  in  low  subsonic  flow.  The  effect  is  slightly  reduced 
when  the  flow  speed  is  increased  and  it  drops  drastically  in  the  transonic  domain.  In  the 
low  supersonic  region  the  force  vanishes  and  does  not  reappear  for  higher  speeds. 

The  side-force  phenomenon  is  therefore  primarily  associated  with  subsonic  speeds.  The 
angle  of  attack  for  onset  of  asymmetric  forces  is  rather  insensitive  to  changes  in  Mach 
number. 


1.3  Reynolds  number 

The  Reynolds  number  may  have  a  dominating  influence  on  the  side-force  induced  by 
asymmetric  vortices.  Keener  et  al  (2)  have  performed  experiments  on  a  tangent  ogive 
forebody  at  Reynolds  number  between  0.3-106  and  3.8-106  (based  on  base  diameter). 
The  Mach  number  was  0.25.  The  onset  angle  for  significant  side-forces  was  insensitive  to 
the  Reynolds  number.  This  is  consistent  with  other  observations  (3).  The  direction  and 
magnitude,  however,  are  strongly  influenced  by  changes  in  Reynolds  number.  Coe  et  al 
(4),  on  the  other  hand,  report  a  small  Reynolds  number  sensitivity.  Deffenbaugh  and 
Koemer  (6)  find  it  difficult  to  determine  the  effect  of  the  Reynolds  number  on 
maximum  side-force. 

There  have  also  been  performed  experiments  with  boundary  layer  transition  strips.  Pick 
(7)  finds  a  significant  reduction  in  side-forces.  This  is  inconsistent  with  the  experiments 
of  Keener  et  al  (2)  where  side-forces  indeed  changed  but  were  not  generally  reduced. 

Jorgensen  (3)  used  a  grit  ring  to  trip  the  boundary  layer.  The  effect  of  this  ring  was 
small  in  the  situations  investigated.  The  lateral  force  as  a  function  of  the  angle  of  attack 
changed  when  the  grit  ring  was  used,  but  the  maximum  force  was  nearly  the  same. 

Increased  Reynolds  number  and  introduction  of  transition  strips  both  promote  earlier 
transition  to  turbulence.  The  various  experiments,  referred  to  here,  give  contradictory 
and  confusing  results.  This  is  probably  because  the  observed  asymmetric  forces  are 
affected  by  the  boundary  layer  in  interaction  with  other  flow  characteristics.  What  is 
observed  when  increasing  the  Reynolds  number  or  introducing  transition  strips  is  thus 
the  result  of  a  particular  experimental  situation  rather  than  an  expression  of  more  basic 
properties. 

The  asymmetries  occur  in  both  laminar  and  turbulent  flow,  so  the  transition  is  presum¬ 
ably  not  an  essential  reason  for  the  observed  asymmetry.  Local  inflectional  instability  of 
the  mean  boundary  layer  profiles  may  be  a  contributes  This  will  be  further  investigated. 
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Adjusting  the  Reynolds  number  will  probably  not,  in  general,  be  a  fruitful  way  to 
eliminate  side-forces. 


1.4  Nose-blunting 

According  to  Keener  et  al  (2)  nose-blunting  reduces  the  side-force  on  both  tangent  ogive 
and  conical  forebodies.  Measurements  on  bodies  with  radius  of  curvature  6.35  mm  and 
12.7  mm,  respectively,  were  performed.  For  a  conical  forebody  the  side-force  reduction 
is  quite  independent  of  the  nose  radius.  It  is  not  so  for  the  tangent  ogive.  The  smaller 
radius  of  curvature  is  found  to  be  the  more  effective  in  reducing  side-forces.  This  is  not 
in  agreement  with  Wardlaw  et  al  (8)  who,  after  considering  all  experimental  data  up  to 
‘975,  conclude  that  blunt  bodies  show  a  lower  level  of  asymmetric  forces  than  sharp 
ones,  but  the  degree  of  bluntness  does  not  appear  important. 

As  indicated  previously,  the  side-force  is  probably  due  to  the  interaction  of  various 
effects.  Some  of  these  are  not  controlled  in  the  experiments.  We  should  therefore  not  be 
surprised  by  the  somewhat  peculiar  results  reported  by  Keener  et  al  (2)  for  the  blunted 
cone. 


1.5  Theoretical  contributions 

There  exist  many  semiempirical  theories  for  separated  flow  around  slender  bodies  of 
revolution.  Angelucci’s  (9)  multivortex  model  is  one  of  the  most  successful.  He  uses 
experimental  separation  point  locations  as  input  in  his  incompressible  and  in  viscid 
calculations.  By  using  the  cross-flow  analogy  (see  Chapter  5)  he  converts  a  three- 
dimensional  stationary  problem  to  a  two-dimensional  unsteady  one. 

A  method  developed  by  Marshall  and  Deffenbaugh  (10)  is  for  the  time  being  the  most 
promising  one-  They,  too,  use  the  cross-flow  analogy,  but  make  boundary  layer  calcula¬ 
tions  instead  of  using  empirical  separation  lines.  Their  model  therefore  incorporates 
viscous  effects.  The  outer  flow  is  modelled  as  a  potential  flow  around  an  expanding 
circular  cylinder  superimposed  on  the  flow  of  a  number  of  point  vortices  describing  the 
wake.  The  inner  flow  consists  of  an  upstream  and  a  downstream  boundary  layer.  Where 
these  separate,  vorticity  is  fed  into  the  inviscid  wake.  There  is  one  empirical  parameter 
used  which  reduces  the  wake  vorticity  compared  with  direct  calculation. 

The  most  apparent  weakness  of  their  model  is  its  two-dimensionality.  The  three- 
dimensional  separation  process  is  more  complex  than  separation  in  two  dimensions  (11). 
A  two-dimensional  approximation  is  therefore  dubious.  Marshall  and  Deffenbaugh  (10) 
make  a  not  too  convincing  separation  line  comparison  with  a  finite  difference  solution. 

We  shall  in  the  sequel  try  to  improve  the  existing  models  by  performing  three- 
dimensional  boundary  layer  calculations. 
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2  POTENTIAL  FLOW 

Prandtl’s  boundary  layer  approximation  will  be  used  in  Chapter  4.  As  input  there,  a 
solution  for  the  inviscid  flow  about  the  body  is  needed.  This  chapter  is  devoted  to  that 
problem. 

If  compressibility,  in  addition  to  viscosity,  can  be  neglected,  a  substantial  simplification 
of  the  equations  results.  We  shall  now  establish  the  condition  for  this. 

In  adiabatic  flow,  the  relation  between  changes  in  pressure  (Ap)  and  density  (Ap)  is 

to  -  to  (2.1) 

where  s  symbolizes  entropy.  According  to  Bernoulli’s  equation,  Ap  is  of  the  order  of 
pv1  in  steady  flow.  The  velocity  of  sound  (c)  in  the  fluid  isy'(3p/9p)s'.  Equation  (2.1) 
can  then  be  rewritten  to  yield 

7*  ~  ©‘  *  (2.2) 


The  fluid  can  be  regarded  as  incompressible  if  A  p/p  <  1,  and  this  implies  a  small  fluid 
velocity  compared  with  that  of  the  sound.  For  Mach  0.25  we  have  A  p/p  ~  0.06,  and 
neglecting  compressibility  should  be  well  justified  below  this  speed.  This  must  not  be 
taken  as  a  strict  limit,  however,  because  experience  tells  that  compressibility  effects  are 
often  negligible  up  to  M  =  0.5  (12). 

In  Chapter  1  we  considered  various  observations  made  on  slender  bodies  of  revolution  at 
high  angle  of  attack.  We  shall,  first  of  all,  look  for  the  origin  of  the  often  large 
asymmetric  forces  which  are  observed.  These  forces  are  most  pronounced  at  speeds 
where  compressibility  is  shown  to  be  of  little  influence  and  diminish  when  the 
transonic  region  is  entered.  Thus  it  follows  that  compressibility  cannot  be  the  cause  of 
the  side-forces.  It  seems,  on  the  contrary,  that  compressibility  may  reduce  the  effect. 
Our  aim  is  to  explain  the  side-force  phenomenon  and  we  shall  make  calculations  under 
conditions  where  the  asymmetries  are  observed  to  be  substantial.  Thus  the  assumption 
about  incompressibility  should  not,  under  the  given  ambitions,  be  a  limitation. 

The  inviscid  flow  can  now  be^obtained^  by  solving  a  perturbation  problem  where  the 
flow  velocity  is  given  by  V  =  V„— V<p.  V-  is  the  free-stream  velocity.  The  perturbation 
potential  (<p)  must  satisfy  Laplace  equation  (12) 

V1*  =  0  (2.3a) 

with  boundary  conditions 


V-n 


on  the  surface 


(2.3b) 


and  the  regularity  condition 

Igradpl  -*•  0  at  infinity.  (2.3c) 

n  is  the  surface  normal. 
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Analytical  solutions  are  known  just  for  a  few  special  geometries.  An  almost  classical 
method,  developed  by  Hess  and  Smith  (12)  will  now  be  outlined.  The  procedure  is  to 
convert  the  differential  equation  to  an  integral  equation  which  is  then  solved  numeri¬ 
cally. 

A  well  posed  problem  is  obtained  if  we  represent  the  solution  by  a  single  layer  potential 
(13) 


*(x)  =  ; 

s 


(2.4) 


r=|x— 1\ 

where  o(£)  is  the  source  strength  over  the  surface  S.  It  is  readily  seen  that  this  potential 
satisfies  the  Laplace  equation  outside  the  body. 


V’  *(X)  =  /  0(15 V*(|)  dSj  =  -  4tt/ a(f)  8  (x-f)dS^  =  0,  for  x*S  (2.5) 


The  boundary  condition  is  given  as  the  normal  derivative  of  the  perturbation  potential 
on  the  surface  and  from  equation  (2.4)  we  have 

*£l-  +/»$  fn(i)dS8  (2.6) 


The  first  term  on  the  right-hand  side  arises  from  the  limiting  process  of  approaching  the 
surface.  When  equation(2.6)  is  inserted,  the  boundary  condition  (2-3b)  yields 

2ira(x)  -  /  ^  (i)a$dS{  =  -n-V.  (2.7) 

This  is  a  Fredholm  integral  equation  of  second  kind  with  the  source  density  o(£)  as  the 
unknown  function. 

The  equation  is  linear  and  we  can  split  the  problem  into  axisymmetric  flow  and  cross 
flow.  These  can  be  handled  separately  and  the  total  solution  is  obtained  by  super¬ 
position. 


Figure  2.1  The  surface  elements 


The  numerical  method  is  based  upon 
approximating  the  body  with  surface 
elements  over  which  the  source  den¬ 
sity  is  assumed  constant  or  a  known 
function  of  one  of  the  coordinates. 
We  are  only  interested  in  axisym¬ 
metric  bodies.  A  considerable  simpli¬ 
fication  is  obtained  when  the  surface 
elements  are  chosen  as  shown  in 
Figure  2.1.  Because  of  the  symmetry 
in  the  problem,  this  choice  is  evident 
in  the  axisymmetric  flow  where  the 
source  density  is  assumed  constant 
over  each  element.  For  the  cross 
flow,  on  the  other  hand,  it  is  not 


obvious  that  an  optimal  choice  of  elements  has  been  made.  It  can  be  shown  (12), 
however,  that  the  source  density  is  a  cosine-function  of  the  angle  with  the  free  cross 
flow  direction.  The  source  density  on  each  element  is  then  chosen  as  a  constant 
multiplied  by  cos  9. 

When  we  on  each  element  choose  a  test  point  xj  (i=l,N  where  N  is  the  total  number 
of  elements)  and  calculate  equation(2.7)  for  each  xj,  N  algebraic  equations  in  N 
unknown  source  strengths  are  obtained.  This  system  is  readily  solved  by  iteration  or 
inversion.  The  procedure  is  carried  through  for  both  the  cross  flow  and  the  axisym- 
metric  flow.  The  results  are  finally  added  to  give  the  total  solution  of  our  problem. 

The  method  has  been  used  to  calculate  the  incompressible  potential  flow  around  a 
cone  with  10°  semiapex  angle.  To  introduce  a  characteristic  length  in  the  problem  we 
have  considered  a  conical  forebody  on  an  infinite  cylinder.  Calculations  have  been 
made  for  the  following  angles  of  attack:  0°,  5°,  10°  15°,  20°,  25°  and  30°.  A 
selection  of  the  results  is  reproduced  in  Figure  2.2. 


Figure  2.2  Inviscid  axisymmetric  and  circumferential  ( dashed  lines)  velocity  for 
various  angles  of  attack  (o)  on  certain  generators 
6  =  circumferential  angle  from  the  windward  side 
V0  =  free  stream  velocity 


Calculations  on  the  same  body  with  a  spherical  nose  have  also  been  performed  and 
some  of  the  results  are  found  in  Figure  2.3. 
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3  THE  COORDINATE  SYSTEM 

For  various  reasons  which  will  be  considered  in  connection  with  the  boundary  layer 
calculation  (Chapter  4)  it  is  advantageous  to  use  streamline  coordinates.  When  equi- 
potential  lines  are  used  as  the  other  set  of  coordinates,  we  have  an  orthogonal  coordi¬ 
nate  system  on  the  body  surface.  The  metric  tensor  is  diagonal. 

The  inviscid  flow  velocities  on  a  cone  at  several  angles  of  attack  were  calculated  in 
the  previous  chapter.  These  results  are  now  used  in  streamline  and  potential  calcu¬ 
lations. 


The  equation  for  the  streamlines  on  the  surface  is 
\Jq  rdfl 
De'"sr 


(3.1) 


where  Uq  is  the  circumferential  velocity  and  U„  is  the  velocity  in  the  generator 
direction.  Furthermore,  r  is  the  local  radius  of  the  body,  8  is  the  circumferential 
angle  and  x  is  the  arc  length  along  the  local  generator. 

The  calculations  are  initiated  either  from  a  stagnation  point  or  from  a  curve  ortho¬ 
gonal  to  the  streamlines.  The  potential  is  given  by 


$  =  /  U  ds 
So 


(3.2) 


where  U  is  the  total  inviscid  velocity  on  the  surface  and  s  is  the  arc  length  along  a 
streamline. 

The  integrations  are  carried  out  by  a  Runge-Kutta  method  (14).  The  streamlines  and 
equipotential  lines  thus  obtained  are  used  in  the  boundary  layer  calculations.  Both 
metric  coefficients  and  curvature  parameters  have  to  be  calculated  during  the  inte¬ 
gration. 

When  the  cone  operates  under  an  angle  of  attack,  the  streamlines  will  deviate  from 
the  windward  side  and  coalesce  at  the  leeward  side  (Figure  3.1).  After  a  certain 
distance  of  integration  the  streamlines  have  to  be  reinitiated  to  retain  the  mesh  size 
within  certain  limits. 


Figure  3.1  Streamlines  and  equipotential  lines  on  a  part  of  the  cone  at  15  degrees 
incidence 


4  BOUNDARY  LAYER  CALCULATIONS 


4.1  The  equations 

A  three-dimensional,  incompressible,  laminar  boundary  layer  flow  is  in  orthogonal 
curvilinear  coordinates  given  by  the  following  equations. 

Continuity: 


1  3u  ,  1  3v 

Hi  3{  5tj 


+ 


3w 

~5I 


-  K,u  -  K2v  =  0 


(4.1a) 


{-momentum: 


u  3u  x  v  3u 
h,  W  h2  5 n 


+  w 


3w 

“3? 


Kjuv  +  K,v2 


(4.1b) 


77-momentum: 


u  3v  I 

h,  3{ 


£  3*  +  w  3£  -  K‘uv  +  KjuJ  =  -  pHT  %  + 


(4.1c) 


where  z  and  w  are  scaled  variables  given  by  z  =>/Re,-z’  and  w  =-y/Re'-wr.  The  dimen¬ 
sionless  lengths  ({,77, z)  and  velocities  (u,v,w')  are  derived  from  physical  quantities  by 
division  by  a  reference  length  L  and  a  reference  velocity  U0,  respectively.  The  metric 
coefficients  hi  and  h2  are  functions  of  {  and  77.  The  parameters  Kj  and  K2  are  the 
geodesic  curvatures  of  the  constant  {-curves  and  constant  q-cuives,  respectively. 
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3h2 

w 


K,= 


1  3h, 
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(4.2) 


Inviscid  streamlines  on  the  surface  and  curves  perpendicular  to  them  are  chosen  as 
coordinate  system.  When  the  inviscid  surface-flow  is  potential,  the  second  set  of 
coordinate  curves  is  equipotential  lines. 

At  the  outer  edge  of  the  boundary  layer  the  momentum  equations  reduce  to 


U  3U  1  3p 

pH7  5T--h7  3? 


(4.3b) 


pKjU1  =  - 


X- 

h2  3 77 


(4.3c) 


where  U  is  the  inviscid  surface-velocity. 

The  boundary  layer  equations  can  now  be  solved  in  their  dimensionless  coordinates  {, 
17  and  z,  but  a  second  transformation  is  advantageous. 

Depending  on  the  pressure  gradient,  the  thickness  of  the  boundary  layer  will  change 
in  the  streamwise  direction.  We  know  from  similarity  analysis  (15)  that  a  constant 
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outer  flow  will  produce  a  boundary  layer  growing  with  the  square  root  of  the  stream- 
wise  distance.  When  the  outer  flow  is  retarded,  the  boundary  layer  grows  faster  and 
with  acceleration,  it  may  even  shrink.  We  shall  make  calculations  in  all  these  domains. 
Thus  it  is  not  obvious  that  a  certain  scaling  in  the  vertical  direction  will  reduce  the 
computation  time.  When  numerical  tests  were  performed,  it  proved  economic  to 
introduce  a  new  vertical  coordinate 

f  =  /((, r\)z  ,  f(S,T|)  =  s'1*  (4.4) 


where  s  is  the  streamwise  distance  from  the  attachment  point.  This  scaling  makes  it 
possible  to  work  with  the  same  vertical  grid-spacing  throughout  the  entire  compu¬ 
tation. 

Equations  (4.3)  and  (4.4)  inserted  into  (4.1)  yield 
Continuity: 


(4.5a) 


(-momentum: 


(4.5b) 


^-momentum: 


u  ,dv  ,  f  af  dv  .  v  ,3v  ,  f  a/  3v  „  ,  r  3V 

K7(3?  +  T  3F  7  3^  3?  >  +  wf  3? -kiUV  +  K2u= *  = 


KjU1  +P 


(4.5c) 


4.2  Numerical  method 

When  phenomena  related  to  high  angle-of-attack  aerodynamics  are  to  be  investigated, 
the  calculation  of  boundary  layer  separation  is  of  great  interest.  The  flow-velocity 
perpendicular  to  the  inviscid  streamlines  (secondary  flow)  will  then,  in  general,  not  be 
small,  and  the  flow  cannot  be  satisfactorily  described  by  a  small  perturbation  of  a 
readily  solvable  problem.  For  these  reasons  we  cannot  use  numerical  methods  based 
on  common  simplifications  such  as  small  secondary  flow  or  small  perturbations. 

Hall  (16)  has  developed  a  finite  difference  scheme  which  relies  on  neither  of  these 
simplifications.  We  shall,  like  Geissler  (17),  use  this  method  in  streamline  coordinates. 


4.3  Discretization 

The  computational  molecule  is  shown  in  Figure  4.1.  The  grid  sizes  in  the  three 
directions  are  A(,  At)  and  Af  respectively.  They  are  kept  constant  throughout  the 


computation.  A  suffix  notation  is  adopted  to  denote  position:  u^m  n 
nAf). 


=  u(lA{,  mArj, 


Figure  4.1  The  computational  molecule 

i 
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The  momentum  equations  are  discretizated  at  the  point  (lA£,mAq,nA£).  Assuming  the 
velocities  known  at  ((1— l)A{,mArj),  (lA£,(m+l)Ar})  and  (lA£,(m— l)Atj)  for  all  n,  we 
shall  calculate  the  velocities  at  ((I+l)A£,mA7j). 


The  finite  difference  approximation  is  given  by 
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-3u. 

(3?)Urn.n 
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2(Af)3  ^vl+l,m,n  vl+l,m,a-l  4  vl-l,m,iH-l 

-2v.-l.m,n  +  v1.1>m>n.1]  +  0(A^+Af3) 


The  accuracy  of  the  discretizated  equations  is  of  second  order. 

We  now  end  up  with  non-linear  equations  for  the  unknown  velocities  at  ((1+1) 
A£,mAi?).  The  equations  are  solved  by  iteration.  When  quadratic  terms  in  the  un¬ 
known  quantities  arise,  we  reformulate  them  in  the  following  way 

.  3uJ+l  ifi  /3uJ  i  /3uJ  , .  ,, 

(u5f)  =  u,<5r>  +u*  (5?)  (4-7) 

where  the  upper  index  j  indicates  the  cycle  of  iteration. 

When  the  momentum  equations  are  discretizated  by  the  approximations  (4.6)  and 
(4.7)  we  arrive  at  the  following  linear  equations 


4  ^n^+l.m.n  4  cnul+l,m,n-l 
a  4  ^  nvl+l,m,n  4  c  nvl+l,m,n-l  _  ^  e 


(4.8a) 

(4.8b) 


The  coefficients  are  given  in  Appendix  A,  and  the  quantities  contained  there  are 
either  known  from  previous  calculations  or  boundary  conditions.  When  velocities  at 
station  4  (see  Figure  4.1)  appear  in  the  coefficients,  we  use  results  from  the  previous 
cycle  of  iteration. 

The  discretization  of  the  equation  of  continuity  is  performed  differently.  The  finite 
difference  approximation  is  evaluated  at  (lA£,mArj,(n— l/2)Af)  as  follows 

(£WV2  -Km.=  -W.,n».n-ll/^+0(AfJ) 

^Vm.n.1/2  =  [U‘+I.m.n  '“l-l.m.n  4ul+l,m,n-l  -ul-l,m,n-l  ]/4A$  +  0(  A*3  +  Af3 ) 

^^l.m.n-1/2  =  ~Vl,m-l,n  4 'W.n-1  -vl,n»-l,n-l l/4Ar?  +  °(Arj3  +  Af3) 

(^?)l,m,n-l/2  =  Kl.m.n  ~ ul+l,m,n-l  +  ul-l,m,n  “ ul-l,m,n-l  l/2Af  +  0(Af3  +  A*3 ) 

(^Vm.n-I/2  =  Krn^l.n  “  Vl,n-1  +  vlm-l,n  ~ vl,m-l,n-l  l/2Af  +  0(Af3  +  Af?3 ) 

uItm,n-l/2  “  Kl.m,n  +  ul-l ,m,n  +  uM,m,n-l  4  ul-l,«n,n-l  l/4  4  °(A^  4  ) 

vl,K,n-l/2  =  tvW,m,n  +  vl-l,m,n  +  vM,m,n-l  4  vl-l,m,n-l  1/4  +  0(A*3  +  Af3 )  (4.9) 

The  equation  of  continuity  can  now  be  discretizated  to  give 


wlm,n-l  +Ar(K,v?+Kav'*v-^?  If  4  2^?  ff 
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and  when  the  streamwise  and  secondary  velocity  components  are  known,  the  vertical 
velocity  can  be  computed  successively  from  the  bottom  to  the  top  of  the  boundary 
layer.  By  discretization  and  use  of  an  iteration  technique,  the  boundary  layer  equ¬ 
ations  have  been  converted  to  a  set  of  coupled,  linear  algebraic  equations. 

The  iteration  procedure  is  as  follows.  The  secondary  and  streamwise  velocities  at 
stations  1,  2  and  3  (see  Figure  4.1)  are  assumed  known,  and  we  start  with  an 
extrapolation 

ul+l,m,n  -  ul>m+l,n  ul,m-l,n  —  ul-l,m,n 

vl+l,m,n  —  vl»na+l,n  Vl»m-I,n  —  vl-l,m,n  (4.11) 

Equation  (4.10)  can  now  be  solved  for  the  vertical  velocities.  We  start  with  the 
normal  velocity  on  the  surface  (which  is  normally  zero)  and  make  a  stepwise  calcu¬ 
lation  up  to  the  outer  edge  of  the  boundary  layer. 

The  vertical  velocities  are  now  inserted  in  equation  (4.8b)  where  extrapolated  values 
(equation  4.11)  are  used  in  the  coefficients.  With  known  boundary  velocities  the  set 
of  equations  thus  obtained  is  solved  by  matrix  inversion. 

A  similar  procedure  is  used  to  solve  equation  (4.8a),  but  now  the  calculated  secon¬ 
dary  velocities  are  used  in  the  coefficients.  This  completes  the  first  cycle  of  iteration, 
and  we  start  the  second  one  with  the  just  calculated  velocities  as  input. 

The  results  are  presumably  improved  by  each  cycle,  and  the  iteration  is  continued 
until  the  results  of  two  successive  cycles  differ  by  less  than  a  prescribed  value.  The 
iteration  is  then  terminated  and  we  proceed  to  the  next  station  in  the  boundary  layer 
where  the  described  procedure  is  repeated. 


4.4  The  coordinate  system 

We  shall  now  motivate  our  choice  of  coordinate  system  for  boundary  layer  calcula¬ 
tions.  The  coordinates  can,  in  principle,  be  based  on  the  body’s  geometry  or  on  the 
inviscid  streamlines.  If  the  body  is  relatively  simple,  a  geometry-oriented  system  may 
be  advantageous.  The  metric  and  coefficients  of  curvature  may  then  be  calculated 
analytically.  If  streamline  coordinates  are  chosen,  both  the  coordinate  system  and  the 
geometry  parameters  must  in  general  be  obtained  by  numerical  integration.  This  calcu¬ 
lation  is  a  major  undertaking  and  it  must  be  repeated  for  every  angle  of  attack 
considered.  Yet  it  turns  out  that  these  extra  calculations  pay  off  in  our  cases. 

A  cone  with  a  spherical  nose  will  be  considered.  The  simplest  coordinate  system,  for 
which  the  geometrical  parameters  can  be  found  analytically,  is  shown  in  Figure  4.2. 
With  a  non-zero  angle  of  attack,  the  stagnation  point  S  is  displaced  from  the  tip  T, 
which  is  the  origo  for  the  body  coordinates.  Becase  of  their  mathematical  structure 
boundary  layer  equations  have  to  be  solved  from  the  stagnation  point,  but  this  is 
impossible  in  these  coordinates.  To  overcome  this  problem,  Cebesi  etal  (18)  intro¬ 
duce  local  streamline  coordinates  in  the  neighbourhood  of  the  stagnation  point.  These 
must  be  calculated  separately  for  each  angle  of  attack  and  be  connected  with  the 
body  coordinates  further  downstream. 


The  initialization  causes  no  geometric  complications  in  streamline  coordinates  (Figure 
4.2). 
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Figure  4.2  a)  Body  and  b)  streamline  coordinates 


The  line  of  three-dimensional  boundary  layer  separation  on  a  cone  at  angle  of  attack 
will  nearly  coincide  with  a  generator  (19).  In  body  coordinates,  it  is  natural  to  choose 
generators  as  one  family  of  coordinates.  Figure  4.3a  indicates  a  grid  of  coordinate 
curves  and  the  line  of  separation.  When  in  the  calculative  procedure  a  point  A  beyond 
the  line  of  separation  is  reached,  the  boundary  layer  approximation  has  broken  down. 
In  Hall’s  numerical  scheme,  results  at  A  influence  the  calculations  at  B,  C,  D,  etc.  The 
hatched  area  is  therefore  excluded  from  integration. 


Figure  4.3  a)  The  line  of  separation  in  body-coordinates  and  b)  streamline  coordinates 


Figure  4.3b  gives  a  typical  situation  for  streamline  coordinates.  During  the  process  of 
calculation,  we  reach  a  point  A'  which  is  beyond  separation.  The  boundary  layer 
calculation  breaks  down  at  this  point.  Calculations  at  B'  are  abandoned,  but  point  C' 
can  still  be  calculated.  This  point  is  upstream  of  separation.  Results  at  D  follow 
before  a  new  breakdown  at  E  .  In  this  way  the  line  of  separation  can  be  traced  down 
the  cone. 

The  mathematical  foundations  for  analysis  of  convergence  and  stability  of  numerical 
schemes  are  well-developed  only  for  linear  systems.  The  results  from  linear  theory  are 
used  as  guidelines  in  nonlinear  problems.  The  confirmation  of  the  validity  of  this 
method  is  dependent  on  numerical  experiments. 

As  far  as  stability  is  concerned,  we  shall  not  perform  analysis  of  the  equations.  Only 
Courant’s  stability  condition  (20),  which  is  developed  for  linear  equations,  will  be 
considered.  This  condition  states  that  the  finite  difference  domain  of  influence  must 
at  least  include  the  continuum  domain  of  influence.  When  only  convective  terms  are 
included  it  yields 
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The  maximum  value  of  the  secondary  flow  v  is  usually  minimized  when  streamline 
coordinates  are  chosen.  This  system  thus  allows  larger  steps  in  the  stream  wise  direc¬ 
tion  (A£)  than  any  system  of  body  coordinates. 

Based  on  these  considerations,  the  choice  of  streamline  coordinates  is  evident. 


4.5  Numerical  test 

We  shall  test  the  numerical  code  by  applying  it  to  a  problem  with  known  solution. 
The  code  is  developed  to  deal  with  boundary-layer  flows  in  three  dimensions.  The 
test  example  should  therefore  have  three-dimensional  character. 

When  streamline  coordinates  are  applied  to  a  three-dimensional  problem,  there  must 
be  a  component  of  secondary  flow  within  the  boundary  layer.  In  applying  the  code, 
we  shall  mainly  work  with  problems  where  the  secondary  flow  in  the  lower  part  of 
the  boundary  layer  is  comparable  with  the  streamwise  velocity  there.  The  numerical 
code  imposes,  in  principle,  no  restriction  on  the  magnitude  of  the  secondary  flow.  A 
proper  test  example  should  therefore  produce  substantial  secondary  flow. 

We  consider  incompressible  flow  over  a  flat  plate  with  inviscid  flow  velocity 

U  =  1 

V  =  -x  (4.13) 

where  U  and  V  are  the  flow  components  along  the  x  and  y  axes,  respectively. 

By  inspection,  symmetries  are  found  which  permit  similarity  simplifications.  The 
boundary-layer  problem  is  then  reduced  to  two  coupled  ordinary  differential  equ¬ 
ations.  When  a  Runge-Kutta  routine  is  applied,  the  equations  can  easily  be  solved  by 
a  shooting  technique. 

It  should  be  noted  that  there  are  numerous  examples  where  three-dimensional  boun¬ 
dary  layer  equations  are  reduced  by  similarity  considerations.  Most  of  them  are  physi¬ 
cally  unrealistic  because  the  inviscid  flow  does  not  satisfy  the  Euler  equation.  The 
chosen  test  example  is  one  of  the  few  which  are  dynamically  possible. 

Hall’s  numerical  method  will  be  tested  against  the  similarity  solution.  Hall  (16)  has 
himself  tested  the  method  on  a  similar  problem,  but  his  numerical  solution  is  more 
straightforward  than  in  the  present  case.  That  is  because  he  calculates  in  Cartesian 
coordinates  where  no  geometric  term  appears.  In  addition  to  that  we  have  performed 
a  vertical  scaling  to  keep  the  boundary  layer  inside  the  initial  grid  system. 

A  set  of  streamlines  is  used  as  coordinates.  If  we  want  to  exploit  the  advantages  of  an 
orthogonal  coordinate-system,  we  must  find  a  second  family  of  coordinates  which  is 
orthogonal  to  the  streamlines.  With  potential  outer  flow,  the  equipotential  lines 
would  have  done  the  job.  In  the  test  example  the  outer  flow  is  rotational  and  we 
choose  instead 


i  =  xe'y 
n  =  y  +  x*/2 


(4.14) 
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as  the  new  coordinates. 

The  geometry  of  this  system  will  now  be  calculated. 


ds1  =  dx2  +  dy2  =  h!2  d$2  +  h22  dr?2 
d£  =  e'y  dx  —  xe"y  dy 
dtj  =  dy  4-  xdx 


These  equations  yield 
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From  equations  (4.16)  we  have 
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The  coordinate  system  used  in  the  numencal  integration  is  shown  in  Figure  4.4.  The 
Runge-Kutta  solution  along  the  lines  A  and  B  was  used  as  initial  values.  The  results  at 
point  P  were  compared  with  the  similarity  solution  there. 


Figure  4.4  The  coordinate  system  used  in 
the  test  example 


•  In  the  first  test  run,  the  mesh  system 
in  Figur  4.4  was  used.  We  then  had 
A£  =  0.02  and  Arj  =  0.1.  Twenty-two 
points  were  calculated  in  the  vertical 
direction.  The  tolerance  of  iteration 
was  fixed  at  10*4.  This  means  that 
iteration  was  terminated  when  the 
maximum  difference  between  one 
cycle  of  iteration  and  the  previous 
one  was  less  than  10*4 . 

Figure  4.5  shows  the  errors  at  point 
P  where  the  calculations  have  pro¬ 
ceeded  six  steps  in  the  flow  direc¬ 
tion.  The  velocity  components  in  the 
streamwise  direction  are  within  2% 
of  the  ’’exact”  results.  The  secon¬ 
dary  component  has  a  slightly  larger 
relative  error  with  a  maximum  devi¬ 
ation  of  2.4%. 

A  reduction  of  the  iterative  tolerance 
to  10*5  produced  only  neglectable 
changes.  The  errors  were  approxi- 
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mately  halved  with  43  vertical  points  (Figure  4.5).  For  neither  the  streamwise  nor  the 
secondary  component  did  they  exceed  1.3%. 


Figure  4.5  Results  of  finite  difference  calculations 

a )  Primary  and  secondary  velocity  profiles  at  P  (see  Figure  4.4) 

b)  Error  in  the  primary  profile 

c)  Error  in  the  secondary  profile 


We  then  returned  to  22  vertical  points  and  reduced  the  step-size  in  the  flow  direction 
to  A£  =  0.01.  This  implied  13  steps  to  reach  P,  and  Figure  4.5  gives  the  effect  on  the 
errors.  A  corresponding  refinement  of  the  grid  in  the  17-direction  gave  insignificant 
changes. 

The  test  runs  prove  that  the  numerical  solution  converges  to  the  ’’exact”  one  when 
the  grid  size  is  reduced.  We  cannot,  however,  draw  any  general  conclusion  about  the 
connection  between  numerical  error  and  grid  size.  This  will  depend  on  each  particular 
problem.  We  have  learned  from  applications,  for  example,  that  large  pressure  gradients 
require  small  steps. 


4.6  Initiation  of  the  boundary  layer 

During  the  discretization  of  the  boundary  layer  equations  (section  4.3),  it  turned  out 
that  as  initial  information  we  need  the  solution  along  two  equipotential  lines. 

We  shall,  as  previously  mentioned,  perform  boundary  layer  calculations  for  cones  at 
various  angles  of  attack.  To  initiate  the  numerical  procedure,  a  boundary  layer  solu¬ 
tion  for  the  apex  region  has  to  be  known.  Because  there  is  no  characteristic  length 
for  an  infinite  cone,  this  will  be  equivalent  to  a  solution  for  the  whole  body.  This  is, 
however,  the  problem  we  intend  to  solve  numerically  and  for  which  there  is  no 
known  simplifying  transformation. 

By  linearizing  the  potential  flow  and  neglecting  curvature,  Howarth  (21)  has  reduced 
the  boundary  layer  equations  in  the  vicinity  of  a  stagnation  point  to  a  pair  of 
simultaneous  ordinary  third-order  differential  equations.  By  using  his  method,  the 
boundary  layer  calculations  for  a  blunted  cone  can  be  initiated. 
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Figure  4.6  I:  The  mvicid  surface  velocity  calculated  for  a  cone  with  spherical  nose 
(radius  R) 

II:  The  experimental  invicid  surface  velocity  for  a  sphere  (R  =  radius) 
with  laminar  boundary  layer 

We  make  a  digression  in  order  to  compare  this  flow  with  the  laminar,  incompressible 
flow  around  a  sphere.  Figure  4.6  gives  the  experimental  inviscid  flow  under  the 
influence  of  the  separated  wake.  The  experimental  point  of  separation  confirms  that 
our  calculated  point  is  reasonable. 

The  only  parameter  with  dimension  of  length  in  the  calculation  of  the  flow  is  the 
radius  of  curvature  (R)  at  the  nose.  We  use  this  length  to  reduce  the  boundary  layer 
equations  to  dimensionless  form.  The  typical  thickness  of  the  boundary  layer  at  the 


10'3  (5  and  R  are  measured  in  meters). 


(4.18) 


One  of  the  conditions  for  the  validity  of  the  boundary  layer  approximation  is  that 
the  change  of  pressure  through  the  boundary  layer  can  be  neglected.  This  is  equi¬ 
valent  to  having  a  thin  boundary  layer  compared  with  the  local  radius  of  curvature  of 
the  body  (15) 


(4.19) 


Even  with  1  mm  radius  of  curvature,  relation  (4.19)  is  satisfied. 
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Another  condition  to  be  met  is  that  the  continuum  assumption  must  be  valid.  When 
the  radius  of  curvature  is  reduced,  the  boundary  layer  shrinks  according  to  equation 
(4.18).  The  continuum  equations  cannot  be  expected  to  hold  when  the  boundary 
layer  thickness  is  reduced  to  a  point  where  it  is  comparable  with  the  mean  free  path 
of  the  gas  molecules. 

We  consider  a  radius  of  1  mm.  Typical  thickness  of  the  boundary  layer  is  then 
10*s  m.  From  the  kinetic  theory  of  gases  (22)  we  have  v  ~  £c,  where  v  is  the  kine¬ 
matic  viscosity,  £  is  the  mean  free  path  and  c  is  the  speed  of  sound  in  the  gas.  For  air 
at  standard  conditions  this  yields  £  ~  10'7  m.  This  implies  that  at  1  mm  radius  of 
curvature  the  thickness  of  the  boundary  layer  is  approximately  one  hundred  mean 
free  paths. 

Even  for  very  small  dimensions,  the  separation  must  thus  be  real  and  not  a  result  of 
applying  the  continuum  theory  or  boundary  layer  approximation  outside  their  do¬ 
mains  of  validity. 

When  the  angle  of  attack  exceeds  ten  degrees,  the  boundary  layer  stays  attached  on 
the  windward  side  making  a  boundary  layer  calculation  possible  downstream  along 
the  body.  We  have  then  obtained  a  continuous  line  of  separation  up  to  which  the 
boundary  layer  can  be  calculated  (Figure  4.7). 

It  is  known  from  experiments  (23)  that  downstream  of  the  nose-region  the  boundary 
layer  is  attached  on  the  leeward  side  of  the  separation  line.  Hence,  reattachment  must 
take  place  to  produce  a  separation  bubble  in  the  nose-region.  For  small  and  moderate 
angles  of  attack  this  bubble  circumvents  the  whole  body.  When  the  body  is  pitched 
to  more  than  ten  degrees,  the  separation  bubble  is  limited  to  the  leeside  of  the  nose 
(Figure  4.7).  The  lateral  extent  of  the  bubble  is  seen  to  be  reduced  when  the  angle  of 
attack  is  increased. 

We  have  thus  worked  out  a  method  of  boundary  layer  calculations  when  the  cone  is 
at  high  angle  of  attack,  but  we  still  face  the  initiation  problem  for  angles  of  attack  up 
to  10°. 

For  a  pointed  cone  there  is  no  known  transformation  which  reduces  the  initiation 
problem  to  ordinary  differential  equations.  When  the  pointed  cone  is  at  zero  angle  of 
attack,  however,  there  exists  a  proper  transformation  due  to  Mangier  (24). 

For  small  angles  of  attack  it  seems  reasonable  that  the  Mangier  boundary  layer  may 
work  as  a  starting  point.  For  moderate  angles  of  attack,  however,  this  seems  more 
doubtful. 

We  shall  now  perform  a  numerical  test  of  the  boundary  layer’s  sensitivity  to  pertur¬ 
bations  in  the  initial  profiles.  The  angle  of  attack  is  chosen  to  be  10°.  In  addition  to 
the  Mangier  solution  a  linear  profile  is  used  initially.  The  difference  between  these 
two  is  large,  but  is  quickly  reduced  when  the  calculations  advance  downstream.  After 
ten  steps  the  difference  between  the  two  calculations  is  less  than  one  per  cent  in  both 
the  streamwise  and  secondary  profile,  except  at  the  top  of  the  boundary  layer  where 
the  secondary  flow  velocity  approaches  zero. 

This  indicates  that  the  effect  of  the  initial  profiles  is  reduced  with  the  distance 
downstream.  The  boundary  layer  will  be  more  dominated  by  the  local  pressure  gradi¬ 
ents  and  geometry.  If  the  Mangier  solution  is  chosen  as  the  initial  profile,  the  calcu¬ 
lated  boundary  layer  will  be  inaccurate  at  the  first  nodal  points.  The  numerical 
example  indicates,  however,  that  the  results  improve  drastically  downstream. 
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This  initiation  is  probably  the  best  we  can  do  for  small  and  moderate  angles  of  attack 
when  blunting  produces  a  separation  bubble  circumventing  the  nose.  For  the  pointed 
cone  the  boundary  layer  theory  is  not  valid  at  the  apex  where  the  potential  flow  is 
calculated  to  be  infinite.  With  the  proposed  method  we  can,  in  principle,  start  the 
calculation  arbitrarily  close  to  the  tip  and  reach  what  is  believed  to  be  a  good 
solution  as  close  to  the  tip  as  we  wish. 


4.7  Results 

Equipped  with  the  initiation  methods  proposed  in  section  4.6  we  have  performed 
boundary  layer  calculations  based  on  the  potential  flows  in  Chapter  2. 

There  have  been  carried  out  very  few  detailed  experiments  on  cones  at  non-zero  angle 
of  attack  in  the  velocity  domain  where  our  assumptions  are  valid.  The  most  relevant 
information  is  obtained  from  Rainbird  et  al  (25)  who  have  made  experiments  in  a  water 
tunnel.  They  observed  the  lines  of  separation  in  an  area  on  the  pointed  cone  where 
the  deviation  from  conical  symmetry  is  small. 

We  have  for  high  angles  of  attack  introduced  blunting  of  the  cone  to  make  initiation 
of  the  numerical  procedure  possible.  This  means  that  we  have  to  calculate  down  to 
the  point  where  the  effect  of  the  blunting  is  negligible.  When  the  distance  from  the 
blunted  nose  is  three  nose-radii,  the  potential  flow  deviates  less  than  1%  from  the 
flow  on  a  similar  but  pointed  cone.  We  terminate  the  calculations  and  compare  with 
experiments  when  the  scaled  boundary  layer  changes  insignificantly  with  further  cal¬ 
culations.  The  scaling  must  be  performed  to  compensate  for  the  downstream  growth 
of  the  boundary  layer. 

Rainbird  etal  have  used  pointed  cones  with  semiapex  angle  5  of  7.5°  and  12.5°, 
respectively.  To  make  comparisons  with  our  calculations  (6  =  10°)  we  have  to  inter¬ 
polate  between  results  obtained  for  the  two  bodies.  This  is  facilitated  by  using  rela¬ 
tive  (a/5)  instead  of  absolute  (a)  incidence. 

The  calculations  were  first  carried  out  with  5  degrees  angle  of  attack.  The  boundary 
layer  thickness  increases  from  the  windward  to  the  leeward  generator  but  stays  at¬ 
tached  all  the  way  (Figure  4.8a).  We  further  notice  that  the  secondary  flow  is  always 
in  the  same  direction.  Rainbird  et  al  observe  no  leeside  separation  for  the  same 
relative  incidence. 

When  the  boundary  layers  are  used  to  identify  separation  their  form  and  not  their 
physical  thickness  is  of  interest.  In  Figure  4.8  there  is  thus  no  vertical  scale  indicated. 

The  two-dimensional  separation  is  a  familiar  and  well-understood  phenomenon.  Under 
calculation,  the  separation  is  localized  where  the  skin  friction  vanishes.  Whew  three- 
dimensional  separation  is  considered  the  flow  Held  is  more  complex  and  the  defini¬ 
tion  of  separation  is  still  under  discussion  (11). 

For  the  present  calculation  it  seems  natural  to  adopt  MaskeU’s  definitions  (26).  He 
has  classified  three-dimensional  separation  as  either  ordinary  or  singular.  We  shall 
consider  the  ordinary  separation  which  results  in  a  free  shear  layer.  It  is  characterized 
by  the  property  that  limiting  streamlines  from  the  mainstream  and  the  separated  area 
meet  tangentially  to  form  what  is  defined  as  the  line  of  separation. 


In  the  present  calculation  we  cannot  go  beyond  separation.  The  identification  of  the 
separation  line  must  thus  be  based  on  information  obtainable  at  the  upstream  side. 
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Figure  4.8a  a=  5“ 
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Figure  4.8d  a=  20° 

6  -  1)  44.7°,  2)  65.1°,  3)  78.8°,  4)  89.8°,  5)  99.7°,  6)  109.0°,  7)  117.8° 


Figure  4.8e  a  =  25° 

0  =  1)  31.6°,  2)  55.8°,  3)  68.9°,  4)  79.2°,  5)  88.5°,  6)  97.6°,  7)  106.5°,  8)  115.5° 

Figure  4.8  Primary  and  secondary  boundary  layer  profiles  in  the  domain  of  the  cone 
where  deviation  from  conical  symmetry  is  small 

The  profiles  are  given  for  a  selection  of  circumferential  angles  9,  measured  from  the 
windward  generator.  It  should  be  observed  that  the  origos  for  the  various  curves  are 
displaced  relative  to  each  other. 

According  to  Maskell  (26),  the  limiting  streamlines  will  here  coalesce  to  form  the  line 
of  separation  as  their  envelope  curve.  The  central  parameter  will  be  the  direction  of 
the  limiting  streamlines  rather  than  the  magnitude  of  the  skin  friction  which  is  essen¬ 
tial  in  determining  two-dimensional  separation. 

When  the  body  is  pitched  to  ten  degrees  or  more,  the  boundary  layer  calculations 
break  down  on  the  leeside.  The  numerical  method  turns  unstable.  Geissler  (17)  claims 
that  this  feature  is  characteristic  of  separation  when  the  calculation  is  effected  in 
streamline  coordinates.  We  want,  however,  to  identify  separation  by  considering  limit¬ 
ing  streamlines. 

The  thickness  of  all  the  calculated  boundary  layers  increases  with  the  distance  from 
the  windward  generator  (Figure  4.8).  The  development  of  the  secondary  flow  is  of 
particular  interest.  Its  magnitude  is  seen  to  increase  when  moving  away  from  the 
windward  side.  A  maximum  is  reached  before  a  more  complex  profile  develops.  The 
inner  part  of  the  boundary  layer  changes  direction  to  form  an  S-shaped  profile.  Under 
these  conditons  the  limiting  streamlines  approach  each  other  in  the  way  assumed  in 
the  separation  criteria.  The  exact  location  of  separation  is  found  by  extrapolating  to 
the  point  where  the  limiting  streamlines  coincide  with  a  generator.  The  skin  friction  is 
small  but  non-zero  here. 

When  approaching  this  special  generator,  the  boundary  layer  grows  rapidly,  indicating 
that  the  main  flow  is  turned  away  from  the  surface.  We  have  thus  identified  a  line 
which  satisfies  our  intuitive  conceptions  of  separation  as  well  as  MaskelTs  formal 
definition. 

Figure  4.9  gives  a  comparison  of  observed  and  calculated  position  of  separation.  For  a 
relative  incidence  of  2.4  the  two  positions  coincide  exactly.  At  lower  incidence  the 
separation  is  predicted  to  take  place  too  early.  At  higher  angles  we  have  the  opposite 
trend. 

There  has  been  included  no  upstream  effect  of  the  separated  flow,  but  this  is  seen  to 
be  of  less  importance  here  than  in,  for  example,  two-dimensional  separation  on  a 
sphere  (27). 
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RELATIVE  INCIDENCE 


Figure  4.9  Comparison  of  observed  (25)  and  calculated  position  of  separation  as 
function  of  relative  incidence 
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5  A  MODEL  OF  THE  SEPARATED  FLOW 
5.1  The  method 

The  inviscid  and  incompressible  flow  around  a  cone  at  various  angles  of  attack  was 
calculated  in  Chapter  2.  These  results  were  then  used  in  the  boundary  layer  calcula¬ 
tions  of  Chapter  4.  The  calculations  broke  down  when  the  line  of  separation  was 
reached.  We  now  intend  to  carry  the  calculations  beyond  this  line.  It  seems  further 
reasonable,  having  obtained  a  first  approximation  of  the  separation  line,  to  attempt  to 
improve  on  it  systematically  by  including  the  wake.  We  therefore  need  a  model  of  the 
separated  flow. 

The  boundary  layer  leaving  the  surface  at  the  separation  line  can  be  modelled  by  a 
vortex  sheet  springing  from  the  surface  of  the  body.  When  this  sheet  is  linked  to  a 
wake  model,  the  boundary  layer  in  the  resulting  external  flow  can  be  calculated  and  a 
new  separation  line  found.  This  iteration  procedure  may  be  followed  to  its  eventual 
convergence. 

The  mathematical  problem  is  surmountable  if  slender  body  theory  is  applied.  This  is 
shown  to  be,  under  certain  conditions,  a  consistent  first  order  approximation  to  the 
fundamental  equations  (28).  These  conditions  are  that  the  stream  is  isentropic  (no 
strong  shock),  the  body  is  slender  and  its  angle  of  attack  is  small.  The  theory  is  thus 
based  on  the  assumption  that  the  body  induces  small  perturbations  on  the  flow.  The 
angles  of  attack,  for  which  leeside  separation  occurs,  are  not  in  this  context  small, 
and  the  cone  perturbs  the  flow  significantly.  The  conditions  on  which  the  slender 
body  simplifications  are  based  are  thus  not  met  and  the  application  of  the  method  is 
questionable. 

The  prime  argument  for  the  validity  of  the  slender  body  method  rests  on  its  ultimate 
success.  The  flow  fields  obtained,  usually  by  visualization  techniques,  are  similar  to 
the  results  calculated  by  for  example  Angelucd  (9). 

We  shall  use  a  modified  version  of  Angelucd’s  method. 

The  flow  separates  from  the  body-surface  along  a  three-dimensional  separation  line. 
The  boundary  layer  is  converted  to  a  vortex  sheet  which  rolls  up  into  a  region  of 
distributed  circulation.  Rotation  is  therefore  present  not  only  in  the  boundary  layer 
and  the  wake  downstream  of  the  body,  but  also  in  the  external  field  which  surrounds 
it. 

The  slender  body  theory  used  here  is  called  the  cross-flow  analogy.  It  is  different 
from  the  standard  version  which  can  be  applied  to  attached  flow  only.  The  time  has 
to  be  induded  when  the  flow  separates  on  the  leeside  and  we  intend  to  model  the 
wake.  The  calculations  are  made  on  an  unsteady  two-dimensional  geometry.  This  is  a 
major  simplification  compared  with  the  original  three-dimensional  flow.  It  is  in  fact 
the  best  we  can  do  with  the  available  computer  resources. 

The  essence  of  the  technique  is  to  introduce  new  coordinates.  We  make  the  following 
transformations 


x  =  x  —  U0  cosa-t 

y  =  y 

z  =  z 


(5.1) 
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where  x,  y  and  z  are  coordinates  linked  to  the  body.  U0  is  the  velocity  of  the  body 
and  t  is  the  time.  The  flow  is  now  considered  in  the  plane  x=0,  see  Figure  5.1.  In  this 
plane  the  flow  Held  is  that  induced  by  a  circular  cylinder  of  time-dependent  radius 
(correlated  with  the  thickness  distribution  of  the  body)  in  a  uniform  free  stream  of 
velocity  V  =  U0  sin  cl 

Vorticity  is  fed  into  the  flow  from  the  point  of  separation.  There  have  been  made 
several  attempts  to  model  the  separated  flow  by  a  concentrated  vortex  core  and  a 
feeding  vortex  sheet  (11). 


t 

a)  THE  CROSS-FLOW  PLANE 


Figure  5.1  a)  The  coordinate  transformation  and  b)  the  body  in  the  cross- flow  plane 


In  the  present  study  the  vortex  sheet  is  represented  by  a  chain  of  discrete  vortices 
(29).  Each  of  these  is  introduced  in  the  flow  close  to  the  point  of  separation  and 
their  subsequent  positions  are  computed  step-by-step  from  the  induced  velocities  on 
the  vortices.  This  model  allows  the  vortex  sheet  to  roll  up  freely  without  the  unphysi¬ 
cal  restriction  of  concentrated  vortex  cores. 

The  Kutta-condition  used  is  different  from  that  chosen  by  Angelucd.  To  establish 
this  we  have  to  investigate  the  vortex  sheet  close  to  the  line  of  separation.  Smith  (30) 
has,  by  simple  arguments,  shown  that  the  vortex  sheet  must  leave  the  body  tangen¬ 
tially.  The  flow  is  shown  qualitatively  in  Figure  5.2. 

With  the  sheet  tangential  to  the  wall  we  can  distinguish  between  an  upstream  and  a 
downstream  side.  On  the  downstream  side  of  the  separation  line,  the  inviscid  stream¬ 
lines  are  nearly  parallel  to  the  separation  line.  From  the  upstream  side,  however, 
the  inviscid  streamlines  on  the  wall  cross  the  separation  line  and  continue  on  the 
vortex  sheet.  The  limiting  streamlines  of  the  viscous  flow  deviate  inwards  from  the 
curved  inviscid  streamlines,  approach  the  separation  line  but  do  not  cross  it.  On  the 
downstream  side,  a  nearly  parallel  flow  is  envisaged  in  which  the  skin-friction  lines  on 
the  wall  and  the  inviscid  streamlines  are  essentially  parallel  (only  the  inviscid  lines  are 
shown  in  Figure  5.2). 

This  implies  that  the  inviscid  surface  velocity  in  the  two-dimensional  model  must  be 
discontinuous  at  the  separation  point. 

For  simplicity  of  calculation  the  vortex  sheet  is  replaced  by  a  chain  of  vortices.  Gose 
to  the  point  of  separation,  however,  a  linear  part  of  the  sheet  is  retained.  The 
procedure  of  calculation  is  as  follows. 
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VORTEX  SHEET 


WALL 


Figure  5.2  Schematic  illustration  of  a  three-dimen¬ 
sional  separation  with  inviscid  stream¬ 
lines  (solid  lines)  and  limiting  stream¬ 
lines  (chain  lines) 


On  the  basis  of  local  in¬ 
duced  velocities  at  the  vor¬ 
tices  at  time  t,  their  posi¬ 
tions  at  t+  At  are  calcula¬ 
ted-  In  the  interval  At  more 
vorticity  is  fed  into  the 
wake.  This  additional  vorti¬ 
city  is  represented  by  a 
linear  vortex  sheet  springing 
from  the  point  of  separa¬ 
tion.  The  sheet  is  replaced 
by  a  concentrated  vortex 
before  the  next  step  of  cal¬ 
culation.  This  one  is  given 
the  same  total  vorticity  as 
the  vortex  sheet  and  is  posi¬ 
tioned  at  the  centre  of  the 
sheet. 

The  two-dimensional  Kutta- 
condition  can  now  be  ex¬ 
pressed  as 


lim  Us  >  0  from  the  upstream  side 

(5.2) 

lim  Us  =  0  from  the  downstream  side 

where  Us  is  the  inviscid  two-dimensional  surface  velocity.  From  this  we  see  that  all 
two-dimensional  vorticity  is  fed  into  the  wake  from  the  upstream  boundary  layer. 
The  Kutta-condition  is  met  in  the  calculations  by  adjusting  the  position  of  the  linear 
vortex  sheet.  Its  inner  end  is  located  a  distance  e  radially  outside  the  separation  point 
and  it  is  directed  towards  the  latest  emerged  vortex.  The  length  of  the  vortex  sheet  is 
chosen  to  be  75%  of  the  distance  to  this  vortex.  The  vorticity  is  assumed  uniform 
over  the  sheet. 

The  sheet  will  induce  singular  velocity  at  its  end  points  (31).  This  is  physically 
unrealistic.  We  are  only  interested  in  the  surface  flow  and  shall  see  that  the  model 
will  give  realistic  conditions  there. 


5.2  The  vortex  sheet 

The  complex  velocity  potential  of  a  vortex  sheet  outside  a  circular  body  of  radius  r0 
can  be  expressed  as  (32) 

*(")— (5-s> 

where  we,  to  preserve  the  boundary  condition  of  no  cross  flow  on  the  body  surface, 
have  included  the  image  vortex  sheet.  The  integration  variable  s  is  the  length  along 
the  sheet  and  y(s)  gives  the  circulation  per  unit  length. 

We  assume  that  the  sheet  is  linear  and  of  uniform  circulation.  Thus 
7(s)  =  k 

ds  =  dOjC  '  *  ( 


5. 
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where  4>s  is  the  angle  between  the  sheet’s  direction  and  the  real  axis. 


The  velocity  in  the  stream  is  obtained  by  differentiating  the  complex  velocity  poten¬ 
tial 


w  =  u— iv 
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which  after  a  tedious  but  straightforward  manipulation  (Appendix  B)  gives  the  surface 
velocity  in  explicit  form. 

The  surface  velocity  induced  by  the  straight  uniform  vortex  sheet  and  its  image  sheet 
will  now  be  examined.  We  choose  unit  body  radius  and  sheet  strength.  The  sheet’s 
length  is  0.2  and  e  -  0.001. 

Figure  5.3  gives  the  results  for  three  different  angles  between  the  radius  vector  at  the 
point  of  separation  and  the  direction  of  the  vortex  sheet.  The  position  of  the  vortex 
sheet  will  later  be  adjusted  so  as  to  meet  the  Kutta-condition  which  requires  that  the 
two-dimensional  velocity  is  discontinuous  at  the  separation  point. 

It  may  be  argued  that  the  velocity  curves,  though  steep  at  the  separation  point,  are 
still  continuous.  Our  calculations  are  fundamentally  based  on  the  boundary  layer 
approximation  and  this  one  will  break  down  in  the  separation  region.  If  we  try  to 
calculate  the  flow  there,  we  shall  go  beyond  the  theory’s  regime  of  validity.  Within 
the  limitations  of  the  theory  the  sheet’s  contribution  to  the  surface  velocity  can  be 
regarded  as  discontinuous. 


Figure  5.3  The  two-dimensional  surface  velocity  induced  by  a  straight  uniform  vor¬ 
tex  sheet  at  various  angles  with  the  radius  vector  at  separation 
6  is  the  circumferential  angle  from  the  point  of  separation 
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We  see  in  Figure  5.3  that  the  upstream  contribution  from  the  sheet  decreases  when 
the  sheet  approaches  the  surface.  This  is  consistent  with  what  is  expected.  At  the 
point  of  separation  the  boundary  layer  leaves  the  surface  and  continues  as  a  vortex 
sheet.  When  the  sheet  is  abruptly  turned  away  from  the  surface,  the  upstream  effect 
of  the  separation  will  be  larger  than  with  a  vortex  sheet  close  to  the  body. 

5.3  The  wake 


In  the  slender  body  approximation  the  calculation  of  the  wake  is  reduced  to  solving 
the  Poisson  equation  in  a  plane  which  is  penetrated  by  the  cone,  see  Figure  5.1.  At 
t  =  t!  (x*ij)  the  cone  is  represented  by  a  circle  of  radius  rj  (=Xitg6,  where  5  is 
the  semiapex  angle)  in  the  complex  plane  a  =  £  +  it).  The  circular  cone  is  equivalent 
to  an  expanding  circular  cylinder  and  the  complex  stream  potential  without  vortices 
(uniform  flow  with  source  and  doublet  at  origo)  is  (31) 

$j(or)  =  ty-osa'-8-  Ino  —  iU0  sina  •  (a  —  (5-6) 

Tq  u 


We  have  already  calculated  the  potential  induced  by  the  vortex  sheet  close  to  the 
point  of  separation 
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Each  of  the  discrete  vortices  in  the  wake  contribute  with 


(5.7) 


(5.8) 


where  Tj  is  the  strength  of  a  vortex  at  ay. 

Because  of  the  linearity  of  Poisson’s  equation  the  various  potentials  can  be  added  to 
give  the  total  potential  for  the  flow. 

There  is  made  no  assumption  about  symmetry.  The  vortices  and  the  two  vortex 
sheets  do  not  have  to  be  symmetrical  about  the  imaginary  axis,  and  this  opens 
possibilities  for  investigating  asymmetric  effects. 

The  vortices  are  introduced  as  point  vortices,  but  they  will  not  retain  this  structure. 
The  vorticity  will  be  diffused. 

Convection  and  diffusion  of  vorticity  are  given  by  Helmholtz  equation  (31) 

Du  =(«-V)V+vV2w  (5.9) 

Dt 

w  =  curl  V 


When  we  neglect  convection  in  our  two-dimensional  model,  equation  (5.9)  reduces  to 
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and  we  impose  the  initial  condition 
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The  solution  of  equation  (5.10)  is 
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with  the  induced  velocity  distribution 

u  <5-12) 

The  total  vorticity  is  constant  but  with  time  it  is  spread  over  an  increasing  area.  This 
effect  is  included  for  two  reasons.  The  vortices  in  the  wake  will  interact,  and  if  two 
of  them  should  get  close  together  they  will  introduce  abnormal  velocities  on  each 
other.  When  a  vortex  gets  close  to  the  body’s  surface,  a  local  unphysical  peak  in 
pressure  distribution  will  result.  These  effects  will  be  damped  by  diffusion. 


The  vortex  flux  in  a  two-dimensional  boundary  layer  will  now  be  calculated.  The 
vorticity  is  given  by 


to 
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(5.13) 


but  the  vertical  velocity  v  can  be  neglected  in  the  boundary  layer.  The  total  flux  of 
vorticity  is  then 


®idy  =  U|/2 


(5.14) 


where  5  is  the  thickness  of  the  boundary  layer.  Thus,  it  turns  out  that  in  a  two- 
dimensional  boundary  layer  the  total  vorticity  flux  is  given  by  the  invisdd  surface 
velocity.  The  actual  boundary  layer  profile  is  of  no  consequence.  We  have  previously 
concluded  that  the  streamlines  on  the  downstream  side  of  separation  are  parallel  to 
the  separation  line.  This  implies  that  the  downstream  boundary  layer  does  not  contri¬ 
bute  to  the  two-dimensional  vorticity  fed  into  the  wake.  All  vorticity  originates  from 
the  upstream  side. 


The  first  order  primary  separation  line  was  calculated  in  Chapter  4.  To  the  same  order 
of  accuracy  we  know  the  invisdd  velodty  at  separation  and  by  equation  (5.14)  we 
have  the  flux  of  vorticity  fed  into  the  wake. 


Figure  5.4  The  calculated  development  of  the  wake  downstream  along  the  body 
a  is  the  angle  of  attack 
a)  a  -  15°,  b)  a  -  20°,  c)  a  -  25° 
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The  calculation  of  the  wake  flow  is  initiated  with  a  linear  vortex  sheet  tangential  to 
the  body.  The  length  of  the  sheet  is  U^At.  U$  is  the  two-dimensional  inviscid  surface 
velocity  at  separation  and  At  is  the  tune  step.  The  distance  e  from  the  surface  is 
adjusted  so  as  to  meet  the  Kutta-condition.  The  further  calculation  is  as  described 
previously. 

We  have  made  calculations  of  the  wake  for  the  cone  at  15,  20  and  25  degrees  of 
incidence,  respectively.  Figure  5.4  shows  how  the  wakes  develop  in  time  (distance 
from  apex).  The  cross  flow  (U0-sina)  converts  the  vortices  downstream  but  the 
circulation  of  the  vortices  opposes  this  effect  and  the  vortices  gather  in  two  limited 
areas  at  the  leeside  of  the  body.  This  is  qualitatively  in  agreement  with  flow  visuali¬ 
sations  (3,4).  They  confirm  that  the  vortidty  is  distributed  over  certain  areas  rather 
than  accumulated  in  one  concentrated  vortex  at  each  side,  as  is  often  assumed  in 
simulations  (33,34).  It  is  also  in  correspondence  with  experiments  that  the  vortices 
are  squeezed  down  to  the  surface  when  the  angle  of  incidence  is  reduced.  At  20 
degrees  of  incidence  and  even  more  pronounced  at  25  degrees,  the  typical  feeding 

sheet  is  lost.  This  is  not  in 
agreement  with  observations 
and  we  shall  return  to  this 
point  in  Chapter  7. 

The  calculations  are  run  up  to 
the  point  where  the  time-aver¬ 
age  surface  velocity  is  essential¬ 
ly  stationary.  We  have  then 
reached  approximative  conical 
symmetry  which  is  reported  to 
exist(19). 

The  induced  surface  velocities 
are  reproduced  in  Figure  5.5 
after  a  time-averaging  of  the 
curves.  Even  though  diffusion 
is  included  to  reduce  the 
short-range  effect  of  the  vor¬ 
tices,  they  produce  unphysical 
local  variations  on  the  surface. 
A  time-averaging  after  the 
approximate  conical  symmetry 
is  reached  gives  curves  without 
abnormal  peaks. 


5.4  Secondary  separation 

When  a  slender  rotational  body  is  pitched  to  an  angle  where  leeside  separation  takes 
place,  several  lines  of  separation  are  observed  (19,35).  The  primary  Unc  is  already 
calculated  (Chapter  4)  and  we  shall  now  deal  with  the  secondary  line  of  separation. 

Figure  5.5  gives  the  circumferential  surface  velocities  calculated  with  the  two-dimen¬ 
sional  unsteady  model.  Combined  with  the  axial  flow  they  give  the  total  surface 
velocities.  This  flow  will  no  longer  be  irrotational  and  there  are  no  equipotential  lines 
available  as  coordinates.  We  therefore  have  to  integrate  lines  normal  to  the  streamlines 
to  obtain  the  second  family  of  coordinates. 

When  the  angle  of  attack  cr  is  15°,  the  attachment  line  falls  along  a  generator  with 
the  angle  9  =  17.5°  from  the  leeside  symmetry  line.  From  this  generator  the  inviscid 
streamlines  deviate  in  both  windward  and  leeward  directions. 


—  y/u. 


— v/u.- 


Figure  5.6c  a  =  25° 

9  -  1)  2.24°,  2)  11.1°,  3)  19.7°,  4)  27.9°,  5)  33.6°,  6)  35.9°,  7)  37.56°,  8)  41.5° 

Figure  5.6  Primary  and  secondary  boundary  layer  profiles  on  the  leeside  of  the  cone 
The  profile*  are  given  for  a  selection  of  circumferential  angle*  6,  meaiured  from  the 
leeward  generator.  It  should  be  observed  that  the  origos  for  the  various  curves  are 
displaced  relative  to  each  other. 
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When  the  body  is  pitched  to  higher  angles,  the  line  of  attachment  moves  towards  the 
leeward  generator  and  they  coalesce  at  a  =  20° . 

The  in  viscid  surface  flow  is  now  used  as  input  for  boundary-layer  calculations.  The 
results  for  a=  15°  are  given  in  Figure  5.6a.  The  boundary  layer  between  the  attach¬ 
ment  line  and  the  leeward  generator  is  attached.  At  the  other  side  of  the  attachment 
line  the  boundary  layer  develops  more  dramatically.  Both  the  circumferential  and 
total  velocity  increase  rapidly,  causing  the  boundary  layer  to  shrink.  The  curvature  of 
the  inviscid  streamlines  produces  substantial  secondary  flow  in  the  boundary  layer. 
The  pressure  gradient  turns  unfavourable  at  8  =  21.5°  and  the  boundary  layer  starts 
growing  before  it  eventually  separates  at  9  =  23.5°. 

At  a  =  20°  the  boundary  layer  develops  smoothly  from  the  leeside  generator  to  the 
domain  with  rapid  acceleration  {Figure  5.6b).  The  boundary  layer  is  then  compressed 
and  a  substantial  secondary  flow  emerges.  The  pressure  gradient  is  adverse  from 
6  =  19°.  First  it  is  very  weak  but  is  increases  slowly  as  we  move  windward. 

When  calculating  the  primary  separation  we  found  that  separation  occurred  in  connec¬ 
tion  with  reversal  of  the  secondary  flow  direction.  It  is  not  so  for  the  present  flow. 
The  direction  of  secondary  flow  is  reversed  at  0  =  32°.  The  boundary  layer  remains 
attached  but  shows  signs  of  incipient  separation  over  a  substantial  distance  down¬ 
stream  of  this  point.  The  imposed  pressure  gradient  is  close  to  the  limit  of  what  the 
boundary  layer  can  resist  without  separation.  Small  perturbations  in  the  vortex- 
induced  crossflow  may  displace  the  line  of  separation  considerably. 

At  25°  the  circumferential  velocity  increases  with  the  distance  from  the  leeside  gene¬ 
rator  and  reaches  a  maximum  at  6  =  31°.  The  pressure  gradient  is  adverse  from  this 
point  on.  Neither  will  the  flow  separate  here  when  the  secondary  flow’s  direction  is 
reversed  (Figure  5.6c).  The  separation  takes  place  at  0  =  42°. 


We  have  integrated  the  inviscid  and  limiting  streamlines.  The  results  are  reproduced  in 
Figure  5.7. 


Figure  5.7a  a  =  5' 
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Figure  5.7 e  a  =  25° 

Figure  5.7  The  calculated  inviscid  (solid  lines)  and  limiting  streamlines  ( dashed  lines ) 
on  the  unwrapped  cone  surface  at  various  angles  of  attack 


The  calculated  lines  of  secondary  separation  are  not  in  agreement  with  the  experimen¬ 
tal  results  of  Rainbird  et  al  (25).  They  report  the  locations  to  be  nearly  independent 
of  the  angle  of  attack.  For  the  present  cone  the  lines  of  secondary  separation  will 
coincide  with  a  generator  approximately  25°  from  the  leeside  symmetry  line. 

We  shall  return  to  this  point  in  Chapter  7  where  these  first  order  results  are  improved 
in  a  second  cycle  of  iteration. 
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6  STREAMWISE-DIRECTED  VORTICES 


It  may  be  interesting  to  investigate  whether  certain  phenomena  which  are  known  to 
occur  in  boundary  layers  will  contribute  to  the  observed  asymmetric  forces. 

We  have  hitherto  made  all  boundary-layer  calculations  under  the  assumption  that  they 
are  stable  and  laminar,  but  it  is  well  known  that  under  certain  conditions  the 
boundary  layer  will  be  unstable  and  undergo  transition  from  a  laminar  to  a  turbulent 
state.  This  is  known  to  have  a  substantial  effect  on  the  location  of  separation  (27). 

It  seems  sensible  to  classify  boundary-layer  instability  as: 

—  Dynamic  (Taylor-Gortler)  instability 

—  Viscous  (Tollmien-Schlichting)  instability 

—  Inflectional  instability 

Dynamic  instability  is  observed  in  a  boundary  layer  flow  over  concave  surfaces  where 
the  velocity  increases  in  the  direction  towards  the  centre  of  curvature. 

When  viscous  instability  is  active,  energy  is  by  viscosity  transferred  from  the  mean 
flow  to  the  disturbances.  The  instability  of  the  laminar  flow  over  a  flat  plate  (Blasius 
flow)  is  of  this  type. 

Inflectional  instability  is  found  when  the  boundary  layer  profile  contains  a  point  of 
inflection.  In  two-dimensional  stream  this  typical  profile  appears  at  an  adverse  pres¬ 
sure  gradient.  That  often  implies  that  we  are  close  to  separation.  For  three-dimen¬ 
sional  stream,  inflectional  profiles  may  occur  even  in  regions  with  favourable  pressure 
gradients.  This  type  of  instability  is  more  effective  than  the  viscous  type  when  transi¬ 
tion  to  turbulence  is  concerned. 

Boltz  et  al  (36)  have  made  an  experimental  study  of  the  transitional  properties  of  the 
boundary  layer  on  a  swept  wing.  For  certain  configurations  they  found  that  a  slight 
increase  in  Reynold’s  number  resulted  in  a  large  forward  movement  of  the  point  of 
transition.  Flow-visualization  studies  showed  that  this  substantial  change  of  the  flow 
was  immediately  preceded  by  the  presence  of  streamwise-directed  vortices  in  the 
boundary  layer.  This  strongly  indicates  that  the  premature  transition  is  caused  by 
stationary  vortices. 

A  few  salient  features  of  the  vortex  system  may  be  noted.  The  vortex  pattern  re¬ 
mained  fixed  in  position  and  the  spread  of  turbulence  in  the  transition  region  indi¬ 
cates  that  the  direction  of  rotation  is  the  same  for  all  vortices.  There  was  found  only 
one  set  of  vortices  and  not  the  double-layer  structure  which  is  indicated  by  Gregory 
eta/ (37). 

For  the  present  analyses  the  main  conclusion  to  be  drawn  from  this  study  is  that 
when  present,  stationary  vortices  in  the  boundary  layer  are  the  dominating  perturba¬ 
tion.  Other  perturbations  are  then  of  less  significance.  When  no  stationary  vortices  are 
observed,  the  transition  is  delayed  until  other  and  less  effective  disturbances  have 
developed  to  turbulence. 

In  the  case  of  viscous  instability,  the  effects  of  Reynold’s  number  on  transition  are 
intimately  related  to  the  intensity  and  frequency  spectrum  of  the  disturbances  present 
in  the  flow  (38).  This  is  believed  to  be  the  case  also  for  other  kinds  of  instability. 

Irregularities  on  the  surface  of  the  body  will  generate  small  perturbations  in  the 
boundary  layer.  The  spectrum  of  disturbances  will  in  general  be  different  for  the  two 
sides  of  the  body  and  streamwise-directed  vortices  may  be  triggered  asymmetrically. 
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The  outcome  may  be  that  transition  to  turbulence  will  take  place  far  earlier  at  one 
side  than  at  the  other  and  in  the  worst  case  we  have  laminar  separation  at  one  side 
and  turbulent  at  the  other.  In  addition,  the  vortices  in  the  boundary  layer  will 
perturb  the  circulation  fed  into  the  wake  at  separation. 

When  ordinary  viscous  instability  is  initiating  turbulence,  the  flow  is  more  stable. 
Small  perturbations  in  the  Reynold’s  number  are  not  known  to  produce  abrupt 
changes  in  the  point  of  transition.  This  is  a  quality  intrinsically  linked  to  the  appear¬ 
ance  of  streamwise-directed  vortices. 


We  are  looking  for  instabilities  in  the  flow  to  explain  the  large  asymmetric  forces 
observed.  One  of  the  candidates,  which  will  be  investigated  further,  is  the  boundary- 
layer  instability  with  streamwise-directed  vortices. 


The  nature  of  boundary-layer  instability  remains  an  exotic  field  in  fluid  mechanics. 
This  is  because  of  the  well  established,  but  not  fully  understood  connection  between 
boundary  layer  instability  and  transition  to  turbulence.  For  any  flow  it  is,  in  prin¬ 
ciple,  possible  to  investigate  stability  properties.  The  mathematical  problem  is,  how¬ 
ever,  so  complex  that  certain  simplifying  assumptions  have  to  be  made.  First,  the 
disturbance  equations  are  linearized  by  regarding  only  small  perturbations.  Thus,  the 
equations  are  valid  only  up  to  the  point  where  neglected  higher  order  terms  become 
significant.  Further,  Stuart  (37)  made  a  major  contribution  to  the  field  by  showing 
that  for  stability  analyses  ’’the  velocity  component  in  the  direction  of  propagation  of 
the  disturbance  may  be  regarded  as  two-dimensional  flow”.  The  problem  is  thus 
reduced  to  the  well-known  Orr-Sommerfeld  equations. 


Inspired  by  Stuart’s  promising  results  we  shall  go  a  step  further  by  neglecting  viscosi¬ 
ty.  A  necessary  and  sufficient  condition  for  instability  is  now  that  the  velocity  profile 
has  a  point  of  inflection.  All  other  profiles  will  be  stable  at  the  inviscid  limit  (39). 

A  fruitful  stability  analysis  under  the  given  assumptions  must  be  made  on  profiles 
with  a  point  of  inflection.  It  follows  from  the  inviscid  equations  that  a  stationary 
vortex  can  develop  in  the  boundary  layer  only  when  the  point  of  inflection  coincides 
with  a  point  of  zero  velocity  (37).  It  is  worth  noting  that  the  calculations  in  Chap¬ 
ter  4  gave  such  profiles  just  prior  to  separation. 


We  choose  the  last  boundary-layer  profile  calculated  before  separation  at  25  degrees 
angle  of  attack  (Figure  6.1a)  to  investigate  the  structure  of  the  streamwise-directed 
vortices  which  may  develop.  The  stability  equation  at  zero  viscosity  then  yields  (37) 
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secondary  flow  V  (Figure  6.1a)  is  well  approxi- 
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where  f(f )  -*•  2^  as  f  -*■  °o  and  o  is  a  numerical  constant 


w  is  the  vertical  (f-direction)  amplitude  and  a  is  the  dimensionless  wave  number  of  a 
periodic  velocity  perturbation. 


-J04-J02  0  J02  .04  V/U0 


Figure  6.1  a)  the  secondary  velocity  profile  and  ■  b)  the  calculated  streamwise-direc- 
ted  vortices 


The  streamlines  in  the  disturbed  boundary  layer  will  now  be  calculated.  If  the  free- 
stream  reference  velocity  is  set  to  100  m/s  and  the  resulting  streamlines  in  the  plane 
of  secondary  flow  are  calculated,  we  get  the  flow  picture  given  in  Figure  6.1b.  This  is 
qualitatively  the  same  vortex  system  as  Stuart  (37)  calculated  for  a  rotating  disk. 
There  are  two  vortex  layers.  The  vortices  close  to  the  surface  rotate  in  one  direction 
and  the  outer  vortices  rotate  in  the  opposite  direction.  In  the  present  case  we  observe 
that  the  inner  vortices  are  extremely  close  to  the  surface.  The  stability  calculations 
have  been  made  without  concerning  viscosity  and  thus  the  velocity  perturbation  does 
not  satisfy  the  viscous  boundary  conditions.  When  the  non-slip  condition  is  imposed 
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at  the  surface  it  seems  likely  that  the  inner  vortex  layer  will  be  wiped  out.  There 
then  remains  a  single-layer  vortex  system  where  all  the  vortices  rotate  in  the  same 
direction.  This  is  in  agreement  with  the  experimental  results  reported  by  Boltz  etal 
(36). 

Stuart  has  shown  in  his  invisdd  analysis  of  stability  that  stationary  vortices  in  the 
boundary  layer  may  appear  only  for  two-dimensional  velocity  profiles  where  the 
velocity  at  the  point  of  inflection  is  zero.  Tobak  (40)  has  carried  out  a  similar 
analysis  without  neglecting  viscous  effects.  He  finds  that  stationary  vortices  may 
appear  for  profiles  differing  from  Stuart's  profile.  The  viscous  corrections  are  small, 
however,  so  the  point  of  inflection  has  to  be  close  to  a  point  of  zero  velocity. 

The  calculated  boundary  layers  in  section  4.7  develop  the  necessary  profiles  just  prior 
to  the  primary  separation.  The  streamwise  profiles  have  also  inflection  points  here 
and  the  secondary  flow  cannot  be  assumed  to  have  such  a  dominating  influence  on 
stability  as  when  it  acts  together  with  only  viscous  destabilizing  effects.  It  is  thus 
unlikely  that  stationary  vortices  will  develop  in  the  calculated  boundary  layers. 

In  Figure  6.1  the  spacing  of  the  vortices  is  of  the  order  of  the  boundary  layer 
thickness.  A  properly  arranged  experiment  should  therefore  reveal  the  vortices,  if 
present.  No  such  observation  is  known  to  the  author. 
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7  ITERATION 

The  first  order  calculation  of  the  flow  around  a  circular  cone  at  non-zero  angle  of 
attack  is  concluded.  The  primary  and  secondary  lines  of  separation  are  localized.  The 
boundary  layer  is  calculated  for  the  whole  body  except  between  the  two  lines  of 
separation. 

We  indicated  in  Chapter  5  that  a  process  of  iteration  might  be  fruitful  to  improve  the 
results.  The  present  results  must  then  be  used  as  input  in  a  second  cycle,  and  so  on. 

A  second  boundary  layer  calculation  where  effects  of  the  separated  wake  are  included 
is  then  performed.  The  Kutta-condition  imposes  a  discontinuity  on  the  invisdd  two- 
dimensional  cross-flow  at  separation.  This  property  is  included  in  the  first  order 
solution  and  the  second  order  boundary  layer  cannot  stay  attached  across  the  discon¬ 
tinuity.  It  is  thus  impossible  to  have  the  separation  delayed  from  one  step  of  iteration 
to  the  next. 

The  calculated  position  of  primary  separation  was  compared  with  experiments  in 
section  4.7.  At  angles  of  attack  below  24°  the  separation  was  predicted  prematurely, 
and  further  calculations  can  only  worsen  this.  For  higher  angles  of  attack  the  itera¬ 
tion  may  improve  the  first  order  results. 

Convergence  in  the  iteration  process  cannot  be  expected  and  we  terminate  the  calcu¬ 
lation  with  the  first  order  results. 

During  the  calculation  of  the  wake,  two-dimensional  boundary  layer  vortidty  was  fed 
into  the  wake  from  the  point  of  separation.  We  have  argued  that  the  downstream 
boundary  layer  does  not  contribute  to  the  vortidty  in  the  wake.  Only  the  upstream 
boundary  layer  has  to  be  considered. 

In  the  first  order  calculation  we  assumed  that  the  invisdd  surface  veloaty  upstream 
of  separation  was  not  influenced  by  the  separated  wake.  Figure  7.1  gives  the  calcu¬ 
lated  first  order  cross-flow  surface  veloaty  when  the  effect  of  the  wake  is  induded. 
The  wake  has  reduced  the  upstream  surface  veloaty  with  approximately  30%  at 
separation.  This  implies  that  the  vortidty  flux  into  the  wake  is  reduced  with  about 
50%  compared  with  the  first  order  solution. 


Figure  7. 1  The  two-dimensional  cross  flow  calculated  for  the  upper  half  of  the  cone 
"  '■  potential  flow 

— — —  — -  first  order  rorticiry  included 

second  order  vortidty  included 
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We  consider  an  angle  of  incidence 
of  25°.  When  the  wake  is  calculated 
with  the  full  first  order  vorticity, 
the  vorticies  gather  in  two  sym¬ 
metric  domains  at  the  leeside  of  the 
body  (Figure  5.4c).  The  typical 
structure  with  a  feeding  sheet  and  a 
limited  area  with  distributed  vorti¬ 
city  does  not  result.  When  the  vorti¬ 
city  is  halved,  however,  the  obser¬ 
ved  structure  of  the  wake  is  ob¬ 
tained  (Figure  7.2). 

When  the  two-dimensional  surface 
velocity  is  calculated  with  reduced 
vorticity  we  get  what  is  reproduced 
in  Figure  7.1.  This  is  quite  different 
from  the  in  viscid  surface  velocity  which  the  leeside  boundary-layer  calculations  in 
Chapter  4  were  based  on. 

Leeside  boundary -layer  calculations  have  been  repeated  with  the  new  outer  flow.  For 
angle  of  attack  25°  and  20°,  separation  will  take  place  at  29°  and  24°,  respectively, 
from  the  leeside  generator.  This  is  rather  close  to  the  observed  secondary  separation 
at  25°  for  both  angles  of  attack  considered  (25)  and  a  major  improvement  compared 
with  the  first  order  results  (section  4.7). 

It  is  now  reasonable  to  try  to  pursue  the  technique  of  Chapter  5  and  introduce 
discrete  vortices  at  both  the  primary  and  secondary  line  of  separation.  The  idea  is  to 
let  the  vortices  develop  a  wake  where  circulation  from  the  secondary  separation  is 
included. 

Secondary  vortices  are  not  introduced  in  the  wake  before  the  primary  vortices  have 
modified  the  wake  to  form  a  basis  for  secondary  separation.  The  outcome  is  that 
primary  and  secondary  vortices  (with  opposite  circulation)  merge  to  give  an  unphysi¬ 
cal  flow  picture. 

This  implies  that  the  technique  with  discrete  vortices  cannot  be  used  in  modelling  a 
wake  which  is  fed  with  vorticity  from  two  lines  of  separation. 

In  Chapter  8  we  shall  investigate  the  stability  of  the  wake  when  the  vorticity  is 
reduced  with  50%  for  all  angles  of  attack  considered. 


! , 
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Figure  7.2  The  calculated  wake  with  second 
order  vorticity 
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8  STABILITY  OF  THE  WAKE 

The  primary  aim  of  this  thesis  is  to  disclose  the  cause  of  the  observed  asymmetric 
forces  discussed  in  Chapter  1.  Certain  perturbations  have  therefore  been  introduced  in 
the  flow  and  their  downstream  development  calculated. 

In  Chapter  6  boundary  layer  instability  was  ruled  out  as  the  main  contributer  to  the 
asymmetric  forces.  Now  we  shall  look  into  the  wake  for  instabilities. 

We  neglect  the  secondary  separation  and  locate  the  primary  lines  asymmetrically.  In 
addition  to  the  geometric  asymmetry  this  also  implies  an  asymmetric  strength  of 
vortidty.  The  boundary-layer  vortex  flux  fed  into  the  wake  decreases  with  delayed 
separation. 

Based  on  experimental  indications  it  is  assumed  that  imperfections  at  the  body’s  tip 
initiate  the  asymmetries.  Thus  only  the  first  ten  point  vortices  are  introduced  asym¬ 
metrically.  We  start  with  the  first  left-hand  vortex  displaced  10  degrees  in  the  leeside 
direction.  Over  the  next  9  vortices  this  displacement  is  reduced  linearly  to  zero.  The 
calculation  of  the  downstream  development  of  the  wake  is  performed  for  angles  of 
incidence  ranging  from  15  to  60  degrees. 

At  15  degrees  angle  of  attack  the  wake  is  stable.  The  nose  perturbation  is  damped 
and  the  wake  develops  with  only  a  slight  deviation  from  symmetry.  This  is  in  agree¬ 
ment  with  the  observations  at  this  angle  (Figure  8.2). 

When  the  body  is  pitched  to  20  degrees  (see  Figure  8.1),  a  small  asymmetry  in  the 
wake  emerges.  The  wake  seems  to  be  so  stable  that  development  of  substantial  asym¬ 
metry  further  downstream  on  a  long  body  is  not  expected.  This  is  confirmed  by  the 
experimental  fact  that  asymmetric  vortex  development  first  appears  when  the  angle  of 
attack  exceeds  twice  the  semiapex  angle  (Figure  8.2). 

In  agreement  with  Figure  8.2,  the  calculated  wake  (Figure  8.1)  is  clearly  asymmetric 
for  angles  above  20  degrees.  The  effect  becomes  more  pronounced  as  the  body  is 
pitched  to  higher  angles.  This  is  again  confirmed  by  experiments  (Figure  8.2).  At  50 
degrees  angle  of  attack  the  vortidty  from  the  two  sides  mixes  considerably.  There  are 
still,  however,  two  distinct  areas  on  the  leeside  with  surplus  of  positive  and  negative 
circulation,  respectively.  They  are  positioned  asymmetrically  and  will  impose  side 
forces  on  the  body.  When  the  angle  is  increased  to  60°,  the  feeding  sheets  are  still 
intact,  but  the  two  domains  with  concentrated  vortidty  are  wiped  out.  The  wake  is 
brought  to  disorder  and  strong  asymmetric  forces  cannot  be  expected.  This  is  again  in 
agreement  with  the  experimental  results  by  Keener  et  al  (Figure  8.2).  They  find  that 
the  time-average  side-force  diminishes  for  all  but  one  Reynolds  number  when  the 
angle  of  attack  exceeds  approximately  60° . 

The  outcome  of  this  stability  analysis  of  the  wake  is  quite  remarkable.  It  not  only 
gives  the  qualitative  development  of  the  wake  as  the  body  is  pitched  through  the 
interesting  range  of  angles.  The  spedfic  angles  where  asymmetries  first  appear  and 
later  diminish  are  quite  well  predicted. 

It  further  confirms  that  small  perturbations  at  the  tip  are  capable  of  triggering  the 
observed  lateral  forces.  The  model  gives,  as  previously  mentioned,  no  possibility  to 
incorporate  the  feedback  of  the  wake  on  the  boundary  layers.  The  state  of  the  wake 
will  influence  the  location  of  both  primary  separation  and  in  particular  the  secondary 
separation.  Thus  the  calculations  cannot  be  expected  to  give  an  altogether  correct 
picture  of  the  flow. 


a)  a-  20 


b)  a-  25 


Figure  8.1  The  calculated  wake  for  angles  of  attack  ranging  from  20  to  60  degrees 
when  a  stnaU  perturbation  is  introduced  near  the  apex 
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Figure  8.2  Experimental  side-force  coefficient  for  various  Reynolds  numbers  as  a 
function  of  angle  of  attack  (2),  M  =  0.25 

CY  ”  side-force/qS,  q  =  f re  e-stream  dynamic  pressure,  S  *  base  area  of  the  cone 


Asymmetric  forces  are  observed  for  both  laminar  and  turbulent  boundary  layers  (3). 
This  must  be  expected  when  the  wake  instability  is  the  main  contributer  to  asym¬ 
metric  forces. 
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APPENDIX  A 

THE  COEFFICIENTS  IN  THE  LINEARIZED  MOMENTUM  EQUATIONS 

ak  =  wl,m  ,n/4-^rA?_0-5/S(AT?)J-(n-1)(UW,m.n+ul.i,m/i)/16s 
-(n-1)g^(vi+iim>n  +  vI-l,m,n)/8h2  s 

bn  =  Kl.m.n  4  ul-l,m,n)/4hi  A$  4  (vl,m+l,n  “  v\^l,n )/4h2  Af» 

-Kl(u1+1,m,n  +  u1.limtn)/4+l/s(Af)* 

~  (n— ^)(vl+l,m,»'l  —  vH,m,n-l  4  vl-l,m,n+l  ~  vl-l,n»,n-l^8b*  s 

<  =  -wl4n,n/Vs'Af  -0-5/s(Af)J  +  (n-l)(uWmn  +u,.l4n  n)/16s 

4  g|  (n-l)(vl+l,m,n  4  vl-l,m,n)/8h*  s 

_  vl-l,m,n(ul+l,m,n  4  ul-l,m,n)/4bl  A^  ~  vI-l,m,n(vl,m+l,n  —vl,m-l,n)/4bj  Al^ 

~  wl,m,n(vl-l,m,n+l  —  vl-l,m,n-l  )/^\/*"Af  4  vl-l,m,n^ul+l,m,n  4  ul-l,m,n)/4 

-  (u1+1  m  n  +  Ui.iim>n)J  /4  +  K2  (u,+l4n  n  +  uU4n  N )2  /4 

4(vl-l,m,n+l  —  ^vl-l,m,n  4  vl-l,m,n-l  )/2(AJ)  s 

4  (n-1KuW,m.n  4  ul-l1m,n)(vl-l,m,n+l  "  vl-l,m.n-l)/16s 

4  5^(n- 1)vl-I,ni,n(vl-l,n,>n+l  -vl-l,m,»-l)/8h*s 
ds 

~  5rj  (n— 1 ) Vj+j >ln >n ( vi+i ,jn,n+l  —  Vl+I ,m ,n-l 8bj  s 

an  =  wl,m,n/4  V^Af  “  0.5/s(Atf  -  (n-l)(uw  m  n+  nUmA)/l6s 
~  35}  ^n_1)  ^vl+l,m.n  4  vl-l,m,n)/8hJ  s 

bn  =  «M^V2h«A^K2(Vl(m^4V1.1,n,n)/44l/s(Ar)1 

(n— ^^ul*l,m,n*l  —  ul*l,m,n-l  4  uM,m,n+l  — ul-l,m,n-l  )A8s 

cn  =  ~  Wl,m,n/V»A?  -0-5/s(Af)J  +  (n-l)(u,+lmill  +  u,.1>m>n)/16s 
4  gjj  (n-l)(v1+limin  +  v,.lm>n)/8h2  s 
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~  (vl*l,m,n  4 vl-l,m,n)(ul,m+l ,n  ul,m-l,n^^2^^ 

~wl4D4»  (ul-l,m,»*l-ul-l,m,n-l)/4‘N/®'Af+K2  ul-l,m,n(vl*l,m,n+  vl-l,m,n)/4 
-Kj  (v1+1>min+  v1.1,m,n)2/44'  (ul*ljn,N-u?-l,m,N)/4h»  ^ 

+  (ul-l,m,n*l  “  ^ul-l,m,n+  Ui.1(inin)/2s(  Af  )2 
+  (n— l)ul-l,ni,n(ul-l  ,m,n*l  ~ul-l,m,n-l  )/^®s 
~ (n—l)ul+l,ni,n(ul+l,m,n+l— uI*l,m,n-l  )/^s 

^)(vl+l,m,n4’  vl-l,m,n^ul-l,m,n+l— ul-l,m,n-l  )/®^2  s 


APPENDIX  B 


EXPLICIT  REPRESENTATION  OF  VELOCITY  COMPONENTS 

The  following  expression  for  the  complex  velocity  emerged  in  section  5.2 


The  endpoints  of  the  vortex  sheet  are  (Figure  B.l) 

i0  uh 

aA  =  r0  e  *  +  Ae  s  (] 

id 

ae  =(ro  +  e)e  ( 

We  are  interested  only  in  the  velocities  on  the  surface  of  the  body.  Thus 
id 

a  =  r0  e  (] 

To  make  the  calculations  more  surveyable,  we  split  equation  (B.l)  and  start  with 

1»  '  Me'"-  - 

In  [cos0  — cos05  —  ^  cos  ip t  + i(sin0— sm0s  —  j^-sin^)]  — 

In  [cos0  — (1  +^)  cos0s  +  i|sin0  —  ( 1  +  ~j  sin  0^]  (j 

With  the  following  substitutions 

A!  =  cos0  —  cos0s  —  -~cos<Ps 
A  2  =  sin0  —  sin0s  — ^-sin05 
A3  =  cos0  —  (1  +  iL)  cos0s 

A4  =  sin0  —  (l  +  j-)  sin0  (] 


Figure  B.  1  The  geometry  of  the  calculation 
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we  get  as  first  term 

&-e  ^sln(-^~  ^-)  =  -jMicos^s+  sin0s) 
2*  xo-o€  2tr  s  s 


{ iln,^i) + iI"n8,xt)  -  “"8(^)j}' = ■ 1 


which  yields 


u,  =  -p-lnt-^-^ylsin^s  [arctg(^)  -  arctg(^)]cos0s 

4ff  A|+  A4  2ir  Aj  a3 

v,  =  --^-ln[^|j-j|]cos0  -^[arctg(-^)-arctg(^)]sin0s 

4ir  A3+A4  iif  Ai  a3 

The  second  term  is 

ik_(!!AZ^€)e10s  =  iL[Ae-ifl_±el(0s'fls‘0)]  =  u2  -  iv2 
2»V  °  '  2ir  ro  ro 

which  implies 

u2  =  — — [AsinS  +  esin(^s  —  Q.  —  0)] 

2irr0 

v2  = - — — [Acos0  —  ecos($s  —  d&  —  0)] 

2tttq 

For  calculation  of  the  third  term  we  start  with 

oo.—tq  »(0-®s)  2 

e  r0(r0+€)e  s  -r§ 

ln<|cos(0-0s)  +  ^cos(0-0s)  -  1  +  i[sin(0-0s)  +  ^-sin(0-0s)l^- 

ln<T(l  +  —  )cos(0  —  0S)  —  1  +  i(l  +  -r-)sin(0  —  0S)> 

L  ro  *0  J 

With  the  following  substitutions 

Bj  -  cos(0  —  0S)  +  ~  cos(0  — 0S)  —  1 

B2  =  sin(0-0s)  +  ^sin(0-0s) 

B3  =(l  +  f)cos(0-0,)-l 
r0 

B4  =  (l+^)sin(0-fls) 


(B.10) 


(B.l  1) 


(B.12) 
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we  get 


2*  ° 


aae-  rg 


^[icos(^s-2^)  -sin(<^s-20)]|^In(||i||)  +  i[arctg(|i)  -  arctg(|^)]|  = 


which  yields 

UJ  "  “  £ lsin(0s-2fl)  -  i-[arctg(|i-)  -  arctg(|i)]cos(0s-20) 
v3  "  ^  “  arctg(  |^)]  sin(0s-  20 )  -  ^  In  [ ||±||  ]  cos(0s-  20 ) 


u3-iv3 

(B.13) 


(B.14) 


The  velocity  components  in  equation  (B.7)  may  now  be  written 


u  =  Ui  +  u2  +  u3 


v  *  v,  +  v2  +  v3 


where  the  right  sides  are  explicitly  given  by  equations  (B.8),  (B.10)  and  (B.14). 
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